from scipy import stats

# Generate random data
data = stats.norm.rvs(size=1000)

# Compute mean and standard deviation
mean = np.mean(data)
std_dev = np.std(data)

print("Mean:", mean)
print("Standard deviation:", std_dev)